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Editor's  note: 

Dr.  Bade,  wVio  joined  this  project  early  in  June,  1954,  had 
to  leave  us  in  September  pursuant  to  action  of  his  draft  board. 
The  present  report,  while  not  in  as  complete  a shape  as  Dr.  Bade 
would  v/ish,  nevertheless  contains  a number  of  contributions  and 
theorems  of  value  for  hydromagnetic  theory.  (W.  M.  E. ) 
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The  geomagnetic  secular*  variations  may  oe  attributed  to  the 

distortion  of  pre-existing  magnetic  fields  in  the  earth's  core  by 

upwelling  motions  of  the  fluid  toward  the  core -mantle  boundary 

(Elsasser,  1950).  The  purpose  of  the  work  reported  here  is  to 

compute  the  hydromagne  tic  effects  of  such  an  up-velllng.  The  main 

problem  is  to  determine  the  effects  at  the  surface  of  the  core, 

2 

since  the  field  in  the  mantle  is  a solution  of  V B = 0 determined 
by  the  boundary  conditions  at  the  surface  of  the  core  and  a t 
inf ini ty . 

Hydromagne tic  affects  in  a homogeneous  conducting  fluid 
are  described  by  the  equation 


V X (v  X B) 


+ V V^'B, 


(1) 


where  as  usual  v"  = {mo)  It  will  be  assumed  that  v = v(r) 

m ' ^ ^ ^ 

is  a given  stationary  velocity  field j we  are  not  concerned  with 
the  problem  of  the  ponderomotive  reaction  of  the  field  on  the  fluid 
motion.  Equation  (l)  is  in  general  difficult  to  solve  exactly. 

Here  a rather  crude  approximation  is  adopted.  Instead  of  letting 
the  Induction  term  V x (v  x B)  and  the  diffusion  term  y 7 B act 


We  use  the  rationalized  mks  system,  see  E]sasser  (1954) 


Concurrently  on  the  f leld  we  let  them  act  consecutively,  '^hat 
Is,  the  problem  is  to  be  solved  approximately  in  two  steps: 

1.  The  magnetic  viscosity  is  set  equal  to  zero  and  the 
fluid  motion^  is  permitted  to  operate  on  the  field  for  a time  t, 

2.  The  fluid  motion  is  stopped  and  diffusion  is  permitted 
to  act  for  an  equal  time  t on  the  field  which  was  produced  by  the 
fluid  motion  in  step  1. 

The  final  field  resulting  from  these  two  steps  is  taken  as 
an  approximation  to  the  field  which  would  be  produced  by  the 
simultaneous  action  of  induction  and  diffusion  for  a time  t. 

In  view  of  the  unprecise  nature  of  t he  ne  thod  used  for 
solving  equation  (1),  and  since  the  upwelling  should  be  fairly 
well  localized,  it  would  be  pointless  to  complicate  the  problem 
by  working  in  spherical  coordinates.  Instead,  the  surface  of  the 
core  will  be  represented  by  the  infinite  plane  z = 0;  the  core 
fluid  will  be  assumed  to  occupy  the  region  z 0,  the  mantle  the 
region  z ^ 0. 

The  local  rectangular  coordinate  system  is  taken  with  the 
z-axis  upward  and  the  x-axis  pointing  south.  Two  types  of  pre- 
existing fields,  upon  which  the  induction  process  acts,  will  be 
cons idered  o 

1,  A poloidal  dipole  field  will  appear  in  the  local 


system  as  a uniform  field 


^x  “ ~ cosv.p. 


(2a) 


Here  is  the  angle  botween  the  positive  z-axls  and  the  direction 

of  the  vector  is  positive  in  the  northern  hemisphere. 

Assuming  that  the  field  is  due  to  a dipole  of  moment  a located  at 


<•  'jj  > 


the  earth’s  center  and  directed  toward  the  south  pole. 


B = ar"^(l  + 3 coa^e)^/^ 

slnup  = (1  + 3 coa^G )“^'^^aln  0 (3) 

coa '4^  = 2(1  + 3 coa^G  )“^/^co3  0. 

Here  r ia  the  radlua  of  the  core  and  6 la  the  colatitude  from  the 
north  pole  of  the  point  under  cona ideration.  In  the  local 
cylindrical  coordinate  ayatem,  the  dipole  field  beconea 
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bp  - “B  sin  C03  <f  Q 

b^  = B sin  4^  ain  (2b) 

b^  = “B  coa  . 

II*  ^ toroidal  field  in  the  y-direction  would  be  repre- 
sented in  general  by 

2 

b^  = b_  = 0,  = B sin  kz^  e"^  , 

A Z y O 

as  can  be  shown  by  separation  of  variables  in  (l),  with  ^=0  and 
the  boundary  condition  b = 0 on  the  surface  z =0.  If  k = 0, 

A-A  O 

the  solution  reduces  to 

\ = P^o- 

Tills  last  form  is  the  one  which  will  be  used  here.  It  should  be 
a good  enough  approximation  to  the  general  solution  near  the  sur- 
face of  the  bore,  ^t  is  stable  against  diffusion.  In  local 
cylindrical  coordinates,  (4a)  becomes 

^ P^o  P^o  °°3<Po» 

^he  pjoloidal  dlpi'olw  field  (2b)  can  be  derived  frcm  the 
vector  potential 


K»VV<*« 


■'it 

T 


ftp  = B(z^3Lr*s-p  “ p^co34>»  cos  <pQ)sin  cp^ 
= B(ZqS1ii  4^  “ PqCOS  4»  cos  <{>q)cos  (p^ 


a„  = 0, 


and  the  toroidal  field  (4b)  can  be  derived  from 
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(2c) 


= “<p  = 0*  = -Qz^p^  cos  9 


0*^0 


\-*yi  ) 


2.  The  Induction  Problem 
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Parker  (1954)  has  shown  that  the  hydromagne  tic  induction 
equation  (l)  (with  = 0)  has  a formal  integral  which  is  exactly 

analogous  to  Cauchy's  Integral  of  Helmholtz’s  equation  for  the 
vorticity  in  hydrodynamics.  ^ similar  integral  can  easily  be 
derived  for  the  vector  potential.  It  is  simplest  to  Start  from 
the  integral  theorem 


(5a) 


which  inplies  that 


/ A 


dr  = 0, 


dt  ^ ~ (5b) 

where  the  integration  is  along  an  arbitrary  closed  path.  Here 
d/dt  denotes  the  substantial  derivative,  i.e.  the  svirface  or  path 
of  integration  is  assumed  to  move  with  the  fluid.  It  is  convenient 
to  work  in  general  coordinates,  as  this  ensures  that  the  results 
will  be  valid  in  the  curvilinear  coordinates  which  will  be  ^x3ed  in 
the  Induction  calculations „ Equation  (5b)  becomes 


/ (A^dx 


m 


dxl:)  = 0, 


where  are  the  covariant  components  of  the  vector  potential  at 
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the  location  of  a given  fluid  particle  at  time  0.  and  A are 
thoso  of  the  vector  potential  at  the  location  of  the  same  fluid 
particle  at  time  t.  wow 

dx|l  = dx“  — 


Hence 


A_  = a, 

U1  t 


" •-  3x“^  ax'“ 

where  Is  an  arbitrary  function  of  x^,  t.  ^^hs  gradient  term  will 
contribute  nothing  to  B,  and  can  be  omitted s 

^ 5-  <«> 

The  Cauchy-Parke r Integral  for  B can  be  obtained  either  by 
taking  the  curl  of  (6)  or  by  direct  calculation  from  (5a).  The 
latter  method  was  used  by  Parker  In  his  original  derivation..  Here 
the  former  method  will  be  used,  as  It  Is  somewhat  shortert 
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utilizing 

; the  definition  of 

a determinant  In  terms  of  the  € •'  . one  finds 
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The  quax-ititles  are  the  contra  variant  components  of  the  field  at 
the  location  of  a giveix  fluid  particle  at  time  0,  and  the  are 
those  of  the  field  at  the  location  of  the  same  particle  at  time  t, 
Let  be  a system  of  Cartesian  coordinates,  T^en 

8x“  3x" 


^ij  “ Smn  3^1 


and  hence 


:^c.- 


dx“ 


- T/2 

_ g / . 


Now  the  quantity  is  known  in  the  Lagrangian  form  of  hydro- 

dynamics to  have  the  value  unity  in  case  of  incompressible  flow. 

This  quantity  is  Just  J {7b)s 


3x 
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8x 
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3x 


n 


dx 
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Q- 


= r§-t 

V e>o/ 


/2 


ox 


m 


= J, 


Hence  J = 1 for  incompressible  flow. 

It  will  prove  convenient  to  work  in  cylindrical  coordi- 
nates p,  z,  where  ^6)  and  (7)  become 
r 


^Po  ^^o 


1 ^Po 


Po 


U.  V » _ 1.  ^ _ 


1 ““o 

ip  p df  “z 
3p  3p  3z 


(8) 


^ B = i 
P J 


|£_b  +i-|£-b  +1^ 

3p^  P P«  WT  <P  3z_ 


Po  5<Po 


b +i-|£-b„  + |^b 

P P_  B«p^  f 3z^ 


+ b„  + -|2—  b 


(9) 


The  equations  (6)  to  (9)  are  based  on  the  Lagrangian  form 


-.*7. 

I 


of  hydrodynamics,  in  which  the  motion  is  described  in  terms  of  the 

1 j i 

trajectories  x = x t)  of  the  individual  fluid  particles. 

However,  it  is  easier  to  investigate  the  character  of  the  fluid 
motion  as  a whole  when  using  the  2ulerian  method  of  velocity  fields. 
For  this  reason,  we  begin  by  using  the  latter  method  and  siiift  over 
to  the  Lagrangian  viewpoint  only  when  we  have  formulated  the  des- 
cription of  a suitable  fluid  motion.  It  will  be  assumed  that  the 
flow  is  incompressible,  V • ^ = 0.  •‘‘he  formula  for  the  divergence 

in  C77lindrical  coordinates  is 

T 1 6v„ 

’ ' X = •?  Ip  (PV'  ^ ^ 


Hence  a velocity  field  defined  by 

Vp  R(p)Z'(s),  v^  = 0,  v^  = ^ R'(p)Z(z) 


(10) 


is  a case  of  incompressible  flow  exhibiting  rotational  symmetry 
about  the  z-axis.  Here  arconstant  and  R and  Z are  arbitrary  differ- 
entiable functions.  Since  v^  = 0,  the  streamlines  all  lie  in 
planes  <P  = constant,  ‘^he  equation  of  the  streamlines  is  R(p)Z(z) 

= constant.  ^Ince  the  fluid  motion  is  steady,  these  streamlines 
are  also  the  trajectories  of  the  fluid  particles. 

In  order  to  get  a tncdel  which  is  physically  reasonable  and 
which  matches  the  problem,  it  is  necessary  to  impose  some  conditions 
on  the  functjons  R,  Z.  It  will  be  assumed  that 

1.  the  fluid  does  not  ci'oss  the  surface  z = 0«,  hence 

v„  = 0 at  z = 0,  or  Z(0)  = 0. 
z 

2.  on  the  z-axls , there  is  no  radial  motion^  hen6e 

Vp  = 0 at  p = 0,  or  R(p) — *0  as  p — » 0. 

These  two  assijimrtions  are  essential.  One  can  also  impose  the  less 
important  but  still  plausible  assumptions  that 


y '-'.I  ■ 


M ■'■ii’T’  ■'■nj  • 


f 

I 

\ 


r 

? 


,.  p'.; 


I 

■ 

V- 


3r 


p 

i 


-8- 


3.  the  only  radial  motion  is  the  spreading-out  as  the 
field  approaches  the  plane  z = 0>  hence  Vp  = 0 at 
z = “oo  , or  Z*(-oo)  = 0. 

4.  in  any  plane  z = constant,  the  motion  at  large  distances 
from  the  point  p = 0 resembles  that  from  a point  source j 


hence  Vp/^p 


.-1 


at  p “ 00  , or  R(oo  ) = constant. 


It  is  easy  to  find  functions  R,Z  satisfying  these  conditions. 
However,  it  is  simplest  not  to  specialize  these  functions  at  this 
stage,  but  rather  to  carry  the  calculation  through  (approximately) 
in  terms  of  the  symbols  R,  Z themselves.  We  will  lise  power  series 
expansions  in  t,  discarding  terms  of  higher  than  the  second  degree. 
From  (lO)  , in  this  approximation. 


pQ  = p + ^ RZ't  + j^RR’Z'^-RR'ZZ'*-  i R^Z'^Jt^ 


= * 


2 r 1 ( 

- ~ E'Zt  + R'^-RR“  + i RR'  ZZ't' 

P 2p^  L P J 


(11) 


and  correspondingly 


p = p-  - RZ’t  + ^ RR’Z'^-RR'ZZ'*- 
^o  _ 2p“L 


^ R'Zt  + ^fR’S.RR"  + i 
ro  ep-- 


(12) 


RR 


zz't^ 


In  (11),  R,  Z.  and  their  derivatives  are  functions  of  the  final 
coordinates  p,  z of  the  particle;  in  (12),  they  are  functions  ^f 
the  initial  coordinates  p^,  z^.  If  one  plans  to  use  (8)  in  com- 
puting the  final  field  one  calculates  the  derivatives  8Xq~/ux- 
as  functions  of  p,  z directly  from  (11): 

= 1 + £ (R'_  J R)Z't  + ^^fR'^Z‘^+  RR'*Z‘^-R'^ZZ'* 

P P 2d^L 


•M  . a 0.0 


- RR"ZZ  - ^ RR'Z’‘-+  ~ RR‘ZZ  ' + -* 
P P ^ 


iT  o 


(13) 


'ij' 
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8z 


° = § Rz”t  + [rh'Z 'z"-RR’Zz’‘'~  f R^Z'z”']t^ 

P 2P"-  I-  ^ -J 


- § (R*-  -i  R’)Zt  -f  ^[r'R^-Rr"-  i R’^+  I RR'*-  RR 

H P 2p^‘-  P P 


dz. 


''j  ZZ  ’ t‘ 
(13) 


= 1 - •-  R*Z't  + ^^|^R»^-RR*'+  ^ 


i RR'1  (Z  Zz")t^ 


3p„  82^  8<p„  8q>^ 

■^o_  o_  o_  'O 


8<P„ 


= = Sf=l* 

On  the  other  hand,  if  one  plans  to  use  (9)  in  computing  one 
must  calculate  the  derivatives  dx^/Sx^*^  from  (12)  , and  then  elimin- 
ate p^  and  Zq  from  the  resulting  expressions  by  means  of  substi- 
tutions of  the  sort 


i-  = i RZ't 


Po  P p» 


R(Pq)  = R(p) [l  + p R' (p)Z' (z)t| , 


which  give  the  final  result  correct  to  terms  in  t*^.  One  finds 

,2 


1^  = 1 - ^(R'-  R)^’t  + ^ [(R’^-RR*)  (Z'2+ZZ“)  - R^Z’^^Jt^ 


= - § RZ*t  - ^ RR*  (Z'Z*-Zz'‘it^ 


2p‘ 


dz 


= S (R**-  R')Zt  + ^[^Rr"' -R'r"+  i R'' 

^ RR'jzz't^ 


(14) 


- i Rr'*  + p 


||r-  = 1 + X R'z't  + ^-,3  (R*^  + RR”  - i RR*)(Z'^  - ZZ'*)t^ 


2P 


Op  = P z _ d£_  ^ d£_  ^ ^ 
^ cJf^  ^ ^ ' 


^ 0,  ^ = 1. 

o o o 

Substituting  from  (l3)  into  (s)  and  taking  the  curl,  or  sub- 
3tlti:ting  from.  ( 14 ) into  (9)  (with  J = 1 since  the  flow  is 


^ -V  , 
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-10- 


incompreaaible ) , one  finds  in  the  case  of  an  initial  poloidal  field 
(2)  that  the  field  as  modified  by  the  fluid  motion  is  given  by 


= -B  sin  ycosf|l-^  (R’  - -iR)z’t 

+ ^ [(R'^-Rr”)  (Z'^-i-Zz'*)-  R2z'2jt2^  (15a) 

. ri 

+ B cos  RZ'*t  + RR' (Z'z'*-ZZ'")t^j  . 


sin»^  sin<^^l  - “ RZ, 't 


.-.1  I 


(15b) 


RR’ZZ**-  RR’Z’^  + ~ R^Z'^  t^ 


• 2P‘ 


B^  = »B  s 


Kj 


-B  siny.--  (r'*  - -i  R*)Zt 

+ ^[rr’"-R'r'*^^  ^ P ^RR'jzz't^l  (15c) 

-B  cosy|^l+  ^ R’z’t  + ^(R’^+RR”-iRH’ )(Z’^-ZZ")t^^ 


In  the  case  of  an  initial  toroidal  field  (4)  one  finds  for 


the  final  field 


Bp  = p sinf  i'.  - ^ (B.’  - i R)Z't 

+ ~ r (r’2-RR‘’)  (Z’2+ZZ"  ) - R^Z'^lt^  (16a) 

2P  L J . . 

2 ) 

- § R’Zt  + (3R*^  - RR**“  -I  RR')ZZ’t2( 

^ 2P‘^  ^ / 


§ R’Zt  + ~ (3R*^ 
^ 2P‘^ 


B = p cos<o/z  - ^ RZ't  •'r  S^Frr' (ZZ’*-Z'^)+  J R^'^lt^ 

r 'Ip  2p'^L  P 


- ^ R’Zt  + (R’^-RR**  + § RR’)ZZ’t^^ 

P 2P'^  P ) 


(16b) 


P r • •%  V . ' • 
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= P sin  «>(—■  (r"  - ” R')Zt  - ^ (R-R«  - iR»^)Z^t^ 

2 (^P  P p 

+ [rR”»  -R>R”  + i R»^  - ^ RR**  + RR’]zZ't^( 


(16c) 


If  the  diffusion  process  is  alow  in  r.omparlson  with  the 
induction  process  (l.e.  if  is  sufficiently  small  in  the  region 
z<'0)j  then  (15)  and  (16)  already  give  a crude  approximation  to 
the  solution  of  equation  (1).  i>ince  Z(0)  = 0„  it  is  apparent 
that  the  field  (16),  produced  by  action  of  the  fluid  motion  on  the 
initial  toroidal  field  (4)^  is  zero  on  the  boundary  z = 0.  Thus 
the  Induction  term  V x (j^ x ^)  in  equation  (1)  is,  by  itself, 
incapable  of  causing  the  toroidal  field  to  manifest  itself  outside 
of  the  core. 

The  condition  Z(0)  =0  does  not  bring  about  such  a marked 
simplification  in  the  field  (15),  which  is  produced  by  action  of 
the  fluid  motion  on  the  Initial  poloidal  field  (2),  A few  terms 
vanish  in  (15a)  and  (15h)^  and  the  first  two  lines  of  equation 
(15c)  disappear.  In  order  to  gain  a certain  insight  into  the 
meaning  of  these  formulas,  let  us  compute  the  average  values 


^(r)  = (itr^)”^  / / Bjj^(P,(p,0)pdpd(p 


r 2ti 


(IV) 


o o 


of  the  Cartesian  comoonents  B,  of  B over  a circle  of  radius  r in 
the  plane  z = 0,  These  are 


= - B 3 


sin  + B sin  vj/*  » ^ ■^o 

- I"  (R'^  - RR"  - ~ RR*  + ^ R^)S'^'(0)t^  I dp 


(ISa) 


B - 0 
7 


( 18b ) 


T* 

B^  = -Bco3i^-2Bco3'^  • aR  ' Z ’ ( 0 ) t 

2 ■ (18c) 

+ (R'^  + RR“  - -i  RR«  )Z«^(0)t“  j dp. 

w —2 

Thus  + r i'^(r),  where  ^ Is  the  initial  field  (2a)  and 

the  quantities  1.  (r)  are  the  changes  in  the  components  of  b due 
to  the  action  of  the  fluid  motion,  integrated  over  a circle  of 
radius  r about  the  origin  in  the  plane  z = 0. 

In  the  mantle  (z  > O),  the  conductivity  a is  small  and  hence 
we  can  set  = oo  , approximately.  Then  in  this  region  B is  a 

Q 

solution  of  V B = 0,  matching  the  field  of  the  core  at  the  boundary 
z ~ U and  reverting  to  the  uniform  poloidal  field  (2)  at  infinity, 

A field  whose  average  (17)  is  zero  over  every  circle  about  the 
origin  in  the  plane  z = 0 will  vanish  more  rapidly  as  z — » oo  than 
a field  whose  average  is  in  general  non-zero.  In  fact,  it  can  be 
shown*^  that  Ij^(oo)  can  be  taken  as  a rough  measure  of  the  amount 
of  field  which  will  penetrate  to  the  surface  of  the  mantle  and 
produce  an  observable  variation  in  the  field  there; 

I^  = B sin  ^ I aR(oo  )Z>  (0)t  - a^Z’^(0)t^  p“^R'^dp|  (19a) 


1=0 

y 


I^  = -2B  cos  aR  ( 00  ) Z * ( 0 ) t j . 


(19b) 

(19c) 


Here  the  conditions  2 and  4 (pp.  7-8)  on  R(p)  have  been  used  to 
simplify  and  evaluate  some  of  the  integrals. 

According  to  this  model,  the  east  horizontal  component  of 
B at  the  surface  should  not  change  appreciably,  and  to  linear  terms 


See  Appendix  A, 
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in  t the  percentage  change  In  the  vertical  component  should  be 
twice  that  In  the  horizontal  component  and  should  have  opposite 
sign.  It  should  perhaps  be  emphasized  that  the  above  conclusions 
have  only  been  demonstrated  in  the  approximation  In  which  cubic 
and  higher  terms  in  t are  neglected,  and  in  which  the  effects  of 
diffusion  are  Ignored. 

3.  The  diffusion  Problem. 

The  problem  of  computing  the  effects  of  diffusion  of  a 
magnetic  field  across  a surface  of  discontinuity  in  o is  in  general 
a formidable  exercise  in  analysis.  In  this  section  we  will  give 
this  problem  a precise  formulation  and  attempt  to  draw  some  con- 
clusions about  the  diffusion  of  the  fields  (15)  and  (15). 

When  V = 0 in  (1),  that  equation  reduces  to  the  diffusion 

AAA 

equation 

P V^B.  (20) 

at  ~~ 

This  eq\iatlon  is  valid  throughout  any  region  where  a = constant. 

If  a Is  a function  of  position,  an  additional  term  a ^Va  x E 
appears  on  the  rlght^,  and  if  a is  discontinuous  at  any  point, 
equation  (l)  is  not  valid  at  that  point. 

If  a is  constant  throughout  space,  the  magnetic  diffusion 
problem  can  be  solved  by  a method  exactly  analogous  to  that  used 
in  the  corresponding  scalar  diffusion  problem  (e.g.  that  of 
conduction  of  heat  in  an  infinite  homogeneous  medium).  In  this 
case,  the  problem  can  be  split  into  three  independent  diffusion 
problems  for  the  Cartesian  components  of  Using  the  method 


- 


. 

I 
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of  Fourier  transforms  > one  finds 


J ^ 00  " ^5  H 

B (r,t)  = (4'itvt)“^' - / B^(r>  ,0)  exp  - d^r ’ . (21) 

-00  ^ — I 4ut  J 


B r (r  - r • 

exp  - 

^ 4i/t 


(r  - r 
4vt 


t ^2 

J 

>t 


and  from  V • ^(r^',0)  = 0,  integrating  by  parts,  one  can  show  that 

the  vector  field  defined  by  (21)  is  solenoidal,  V • _^(r,t)  = 0» 

The  same  result  follows  directly  from  BB/dt  = - x (V  x B). 

^ ^ ' in 

The  problem  under  consideration  here  is  more  intricate,  as 
a is  discontinuous  at  the  surface  z = 0: 


0=0^  for  z < 0 

0 = 0 for  z 0 


Consequently; 


^m  = ^ = (A<o/ 


V =00 
*^m 


for  z < 0 
fOx*  z > 0 


o 

Hence  in  the  region  z 0 , ^satisfies  the  equation  ^ B^  = 0.  At 
z = 0,  equation  (20)  is  not  valid.  According  to  the  analogy  with 
scalar  diffusion  theory,  one  would  seek  to  determine  six  functions 

^xl * ^y>’  ^zl’  ®x2 ’ ^y2’  ^z2  position  and  time,  such  that 
B^(r,0)  is  the  initial  field  in  z <•  0,  B„(r,0)  is  the  initial  field 
in  z )>  0,  the  vector  _B^  , t ) satisfies  equation  (20)  with 
in  the  region  z < 0,  and  the  vector  _Bg(^,t)  satisfies  V _Bg  = 0 
in  the  region  z > 0.  Since  the  diffusion  equation  does  not  hold 
at  z = 0,  the.  behavior  of  the  functions  on  the  boundary 


See  Sneddon  (1051),  chapter  5. 
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must  be  determined  from  boundary  conditions.  There  should  be  six 
of  these,  as  there  are  six  functions  to  be  determined. 

In  order  to  find  these  boundary  conditions , one  may  first 
note  that  the  behavior  of  the  field  at  the  boundary  is  certainly 
determined  by  the  f xjindamental  equation3  of  electromagnetic  theory 
in  a medium  at  rest: 

V X ^ = - 8B/dt,  eV  ‘ ^ “ P (22a) 

V X B = dE/dt  +/<-J,  V • B = 0 (22b) 

J = oE  (22c) 

^ * i.  = 0.  (22d) 

In  diffusion  theory,  the  boundary  conditions  are  relations  which 
must  be  satisfied  for  t>0,  but  not  necessarily  for  t = 0,  »Vhen 
t > 0,  E and  B and  their  time  derivatives  are  certainly  finite. 
Hence  it  follows  from  (22b)  that  at  z = 0,  Iho  vector  ^ is  con- 
tinuous s ~ JS-2’  the  first  equation  of  (22a),  it  follows 

tliat  F and  E_  are  continuous  at  z = 0.  From  the  second  equation 
(22a),  E *•  E , = "c/e  ? where  x Is  the  surface  density  of  charge 
on  z = 0,  According  to  (22d),  J^g  - = - 8Tr/3t.  Combining 

these  last  two  relations  with  the  first  equation  of  (22b),  one 
finds  tb^  t (V  x B)„  is  continuous  at  z = 0. 

In  the  approximation  usually  used  in  hydromagnetlc  theory, 
the  displacement  current  e ^/3t  is  neglected  in  comparison  with 
the  conduction  current  In  this  same  approxJna  tion,  the 

quantity  dT/dt  = e(8E^g/6t  - 8E^^/dt)  must  also  be  neglected. 

But  then  >x  ^^7  t>e  regarded  as  constant  in  time,  and  the  corres- 
ponding longitudinal  part  of  ^ may  be  split  off  and  discarded. 

The  justification  for  this  procedure  is  that  the  longitudinal 


-16- 


part  of  ^ will  be  e slowly  varying  electric  field,  and  hence  will 
have  little  effect  on  the  rragnetlc  field  B.  In  this  approximation, 

A/V\ 

J Is  continuous  across  the  boundary.  But  In  the  region  z > 0. 

J = 0 because  o = 0:  hence  J„  = E„  = 0 at  z = 0. 

The  preceding  discussion  has  established  the  six  boundary 
c ondltlons 

®xl  ~ 

^xl  “ ^x2’ 

In  order  to  employ  these  conditions  In  the  solution  of  the  mag- 
netic diffusion  problem  In  the  form  outlined  above.  It  would  be 
necessary  to  express  them  In  terms  of  derivatives 

of  these  vectors.  However,  this  is  not  possible  for  the  conditions 
on  E and  E , because  (22b)  and  (22c)  give  the  only  direct  relation 

y 

between  E and  the  derivatives  of  B,  and  this  relation  breaks  down 
in  the  region  z > 0,  where  o = 0. 

Here  appears  the  essential  complicating  distinction  between 
magnetic  and  scalar  diffusion  problems.  The  boundary  conditions 
of  the  magnetic  problem  may  each  invol’-'e  several  components  of 
Bi,  as  in  a condition  on  V x as  a consequence,  the  problem 
cannot  generally  be  split  Into  three  Independent  problems  for  the 
Cartesian  components.  Moreover,  the  boundary  conditions  are 
conditions  on  the  electromagnetic  field  which  cannot  be  reduced 
to  conditions  on  the  magnetic  field  alone. 

This  last  difficulty,  however,  can  be  circumvented  In  a 
simple  wayj  If  the  problem  is  formulated  In  terms  of  the  po- 
tentials rathep  than  the  field  intensities,  one  can  express  all 


B,=B„,  B^=B„ 

yl  y2’  zl  z2 


(23) 
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Of  the  boundary  conditions  in  terms  of  the  same  dependent  variables 
that  appear  in  the  diffusion  equation.  As  mentioned  a'oove,  if  the 
displacement  current  is  neglected  (as  Is  always  done  in  hydro- 
magnetism),  then  the  quantity  dr/dt  must  also  be  neglected,  and 
the  longitudinal  component  of  _E  due  to  a constant  surface  charge 
T can  be  split  off  and  discarded.  Then  the  field  intensities  can 
be  represented  in  terms  of  the  vector  potential 

B = V X A,  E = -3A/^t  (24) 

where 

V * A = 0.  (25) 

Then  all  the  equations  (22)  are  satisfied  identically  (with 
p = 0 everywhere)  except  the  first  equation  of  (22b),  which  becomes 


^ = V_V^A. 
dt  ni 


(26) 


r- 

k 

i: 


i.W. 


■ 
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Since  ^ is  finite,  it  follows  from  the  first  equation  of  (24)  and 
from  (25)  that 


^xl  ^x2*  ^yl  ~ -^y2»  ^zl  ^ 


z2 


(27a) 


The  former  boundary  conditions  (23)  on  namely  become 

V X Jji  = V X (27b) 

In  view  of  the  second  equation.  (24),  (27a)  implies  that  and 

E are  continuous  at  the  boundary.  The  remaining  condition  (23), 

«/ 

namely  E = 0,  becomes 
z 

— — 0 at  z — 0.  (27c) 

Wo  now  have  exactly  six  Independent  conditions  on  Ac 
Equation  (27a)  gives  three  conditions}  (27b)  gives  only  two  more. 
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since  the  continuity  of  A and  A implies  the  continuity  of 

^ y 

(.7  X A)_t  and  (27c)  gives  the  aiicth  condition. 

The  fact  that  five  of  t^hese  six  conditions  are  statements 
of  continuity  means  that  it  will  be  necessai^  to  solve  simultaneously 
for  the  two  vector  fields  A.  , A^,  where  A.,  (r,0)  and  Ao(i’.O)  are 


given  and 


vJ  7*^^  = 3^/dt  in  z < 0 


7^2  = 0 


in  z > 0 


(28a) 

(28b) 


and  the  boundary  conditions  (27)  are  satisfied  at  z = 0. 

The  solution  of  this  problem  can  be  obtained  by  means  of 
the  following  procedure.  According  to  (27c),  A (x,y,0)  is  not  a 
function  of  t.  In  the  region  z^  0 we  may  set 

A^,(r,t)  = u(r)  + v(r,t),  (29) 


where  v(x,y,0,t)  = 0.  Then 

= 0,  If  = . 

The  solution  of  this  diffusion  problem  for  v is  known  from  the 
theory  of  heat  conduction  ( semi-inf inte  region,  temperature  zero 
on  the  boxxndary) . It  is  p 


(r.,t)  = (47tvt)"^/^  // / v( 


•00  -CD 


-fj_^/4yt  3 

» 0)  e ^ -e  ~ d r « 

> L j ^ > 


where 


Now  X© t 


= (x-xM^  + (y-y>)^  + (z-z«)^ 
fg^  = (x-x«)^  + (y-y*)^  + (z+z»)^. 


J ^ = (x-xO^  + (y-y*) 


2 . 2 

+ z . 
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Then  According  to  equation  (1)  of  Appendix  A,  the  field  u(r) 
in  the  region  z <C  0 is  given  by 

_ OO 

u = - //  ,y » ,0,0)dx»dy'  (3l) 

-00 

and  the  field  A^g  in  the  region  z>0  is  given  by 

^z2  “ 2it  ,y  ' ,C,0)dx*dy  * (32) 

-CO 

The  field  A^  is  now  determined  throughout  space  by  equations 
(29),  (30),  (31),  and  (32). 

The  problem  of  determining  A^  and  A is  more  difficult, 

y 

because  of  the  nature  of  the  boundary  conditions  on  these 

quantities.  According  to  equation  (l)  of  Appendix  A and  the 

boundary  conditions  (27a).  A and  A can  be  represented  in  the 

X y 

region  ?.  > 0 by 


00  - 

^x2  = I?  ,0,t)dx*dy» 

-00 

(33) 

Ay2  = Ay^(x»  ,yS0,t)dx’dy 


The  remaining  boundary  conditions,  which  state  that  3A^/3g  and 
3A  /8z  ire  continuous  at  the  boundary,  must  now  be  used  to  obtain 
the  solution  in  the  region  z < 0,  This  problem  can  be  formulated 
most  simply  in  terms  of  Green's  function,  which  is  here  defined 
as  the  field  G(x,y,z°,  x ’ -y  ’ j z » 5 t)  at  the  point  (x’,y',z*)  at  time 
t due  to  an  Instantaneous  point  source  of  unit  strength  generated 
at  the  point  (x,y,z)  at  time  zero.  It  is  assumed  that  z,z'  <.  0, 
that  the  field  in  the  region  z < 0 at  t = 0 is  zero  everywhere 
except  at  (x*,y*,z*)  and  that  the  boundary  conditions  are 


satisfied  at  z = 0.  Then  one  can  show  that  the  field  at  the  point 

(x,y,z)  at  time  t due  to  an  Initial  field  A (x ’ ,y ’ , z ’ ,0 ) at  time 

zero  is  given  by  | 

oo  o 

Axi(x„y,z,t)  = tT* / A^^(x» ,y » ,z ‘ ,0)G(x,y  ,z;x» ,y » ,z ’ jt )dx’dy »dz » . (34)  | 

— GO  —GO  i 


Thus  if  one  can  determine  Green's  function  G for  the  problem,  (34) 
gives  the  solution  in  the  region  z < 0 for  an  arbitrary  initial 
field  A (x * ,y * , z * ,0 ) . Once  this  field  is  known,  the  field  in 
the  region  z >•  0 is  determined  by  (33). 

If  G(x,y,z,t)  is  any  solution  in  the  region  z <.  0 of  the 
diffusion  equation  and  the  boimdary  conditions  on  A^  and  A 

«/ 

(in  particular,  it  may  be  the  Green's  function),  then  G(x,y,0,t) 
is  the  field  cn  the  boundary.  Substituting  into  (33)  and  applying 
the  bom'xdary  condition  on  dG/oz  , we  obtain 

fdGl  _ 1 ^ rr  z G(x”  ,y"  ,0,t ) , «,  u 

Z-^O  -OD  [(x-x'*)2  + (y-y")^  + z~Y^ 


The  author  has  not  yet  succeeded  in  determining  a Green's  function 
which  aatlsfles  (35)  in  the  general  three-dimensional  case. 

There  is  one  special  case  in  which  the  entire  problem  can 
be  solved  with  ease:  the  case  in  which  A and  A are  functions 

X j 

of  z only.  Then  in  equations  (33.*,  the  quantities  A^^  and  A^^ 
can  be  brought  out  from  under  the  Integral  signs  and  the  integrations 
c an  be  performed  explicitly.  '-‘he  results  are  that  A^g  and  A^g  are 
independent  of  position  in  the  region  z "i?  0.  Thus  3A^/dz  and 
3A^/3z  are  zero  at  z = 0,  and  in  the  region  z < 0 the  diffusion 


■'Carslaw  and  Jaeger  (1947),  chapter  AIII  t 
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problem  for  A and  A is  formally  exactly  analogous  to  the  problem 
of  heat  conducltlon  in  a semi-infinite  iredium  with  the  adiabatic 
boundary  condition.  '^he  solution  for  is  thus 


A^l  = (4Ttvt) 


-1/2 


-(z-z')V4Vt  -(z+z»)^/4vt 


2 


/ A, 


,0) 


dz ' , 


-CO 


and  Ay-  is  given  by  a similar  equation. 

A slightly  different  approach  to  the  general  (three- 

dimensional)  problem  is  obtained  if  it  is  assumed  that  , and 

t'g  are  both  finite,  but  that  2 ^1‘  boundary  conditions 

(27)  remain  exactly  the  same.  AH  components  of  ^ and  v x A ^ are 

continuous  across  ttje  boundary.  From  the  electromagnetic  equations 

(22),  one  finds  that  (V  x B)  must  be  continuous,  and  hence  that 

dA  /dt  must  be  continuous  across  the  boundary.  Bat  A itself 

is  continuous,  so  that  (unless  1 ~ ^3^  conclude  that 

8A  /8t  = 0,  as  before.  Now  the  diffusion  equation  is  satisfied 
z 

both  below  and  above  the  boundary  (but  with  different  values  of  >^  ) , 
so  that  both  f ielda  and  can  be  represented  in  terms  of 
Green's  functions.  It  may  prove  to  be  an  advantage  to  have 
functions  of  the  same  type  on  both  sides  of  the  boundary. 

In  the  composite  medium  with  finite  i'g,  the  boundary 
conditions  can  also  be  stated  in  t erras  of  B rather  than  of  A. 

They  are 


— 1 —2 
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(V  X = (V  X Bg)^  = 0. 
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(The  laat  of  these  holds  in  the  approximation  used  here,  in  which 
is  assumed  to  be  continuous  at  z = 0.  ) In  the  limit  V' g — ^ oo  , 
the  conditions  on  (V  x B)  and  (7  x B)._  in  their  present  form  give 
no  information  about  . However,  these  conditions  can  be  trans- 
formed  in  tbls  case  into  conditions  somewhat  resembling  (35).  In 
the  region  z > 0,  these  conditions  state  that  when  V_  - co  ^ 

V X JBg  = 0,  and  hence  = VV.  Now  since  B^g  = dV/oz  is  given  on 
the  boundary,  V is  determined  throughout  the  region  z > 0 by 


V = - ^//f  -^B  (x»,y»,U)dx»dy'. 


B^.(x,y,0,t)  = B-(x,y,0,t)  = 


are  the  required  conditions.  It  should  be  noted  that  here  the 
various  components  of  ^ are  very  much  mixed  together  in  the 
boundary  conditions.  One  great  advantage  of  the  formulation 
based  on  the  vector  potential  is  that  it  proves  possible  to 
separate  the  problem  into  three  independent  diffusion  problems 
for  A^,  Ay,  Ajj. 

It  is  possible  to  extend  the  preceding  results  of  this 

section  to  the  case  of  a spherical  boundary  between  the  two  media. 

The  boundary  conditions  on  A are  that  A and  V x A are  continuous 

across  the  boundary,  and  that  the  normal  component  6A^/3t  = 0 

2 

at  the  bo’undary  r = a.  The  field  ^^2  “ ® r a can  be 

expressed  in  terms  of  the  field  on  the  surface  by  means  of 
integrals  similar  to  equation  (l)  of  Appendix  Aj  see  the  same 


Jeffreys  (1950),  p.  221. 
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reference  to  Jeffreys.  The  field  A_  can  be  determined  by  a method 

X 

similar  to  that  used  in  the  case  of  A with  a plane  boundary; 
the  solution  of  the  diffusion  problem  in  the  interior  of  the 
sphere  with  v = 0 on  the  surface  is  very  much  more  complicated 
than  the  corresponding  solution  in  a semi-infinite  region,  but 
it  is  given  in  the  literature  . 


Carslaw  and  Jaeger 


(1947),  pp.  210-212 
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Appendix  A:  Solution  of  Laplace's  Equation  in  a 

Seml-lnTinite  Region. 


In  the  case  of  the  semi-infinite  region  z > 0,  the  well- 
known  integral 


f (x 


.y,z)  = Is//[f  " - F • 


dS 


of  Laplace's  equation  can  be  transformed  into  the  expression 

f (x«  .0)dx'dY' 


f (x,y ,z)  = // 


[(X-X')2  ^ (y-y.)2  Z^772  ' 


(Al) 


which  gives  the  values  of  the  harmonic  function  f throughout  z > 0 
in  terms  of  its  values  on  the  plane  z = 0. 

One  application  of  this  formula  is  the  justification  of  the 
use  of  equations  (19)  to  approximate  the  field  at  the  surface  of 
the  mantle  due  to  an  upwelling  at  the  surface  of  the  core.  He 
can  split  the  total  field  ^ into  the  constant  (initial)  field 
and  a variable  field  p produced  by  the  action  of  the  fluid  motion 
on  b.  We  assume  that  the  upwelling  is  sxifficiently  well  localized 
that  the  improper  integrals 


^ oo  27t 

/ / p.  (p,  ,0)pdpd<p  (A2) 

o o 

exist.  (This  will  be  the  case,  for  instance,  if  condition  4, 
p,  8,  on  R(p)  is  satisfied.)  Now  the  variable  field  at  the 
surface  of  the  mantle  is  given  by  (Al)  with  Pj^(x,y,z)  substituted 
for  f(x,yjz)5  here  z is  regarded  as  representing  the  height  of  the 
surface  of  the  mantle  above  that  of  the  core.  A representative 


*A  proof  is  given  in  Jeffreys  (1950),  p.  221. 
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Multiplying  equation  (l)  by  3,  one  obtains 


i 3(B^)/3t  = B • 


t; 


- V„B,-  [v  X (7  ^ B)l.  (Bl) 

L J 


V X {v  X B) 

The  first  term  on  the  right  can  be  transformed  using  vector 
identities  in  the  following  ways 


B 


fv  X (v  X B)1  = B • f(B  • V)v  - (v  • V)B  - B(V  * v)] 

1^  J I ..v%\  — ^ **%•>»  I 

= B.B,  dv,./dx.  - -5  V • V(B^)  - B^(y  ' v). 

i J 1 J CS  /w\ 


Le  t 


, /Bv,  , , 

V --  1 / __i  + _J  ] 

IJ  c'Xj  dx^ 


(B2) 


Then 


B • [V  X (v  X B)]  = “ 5 Z • 

The  second  term  on  the  right  in  (Bl)  can  be  transformed 
as  follows: 

• [Vx  (VxB)]=7  • [b^x  (Vx^)l  - (Vx^)^ 

= ~ V^(B^)  - (V  X B)^  - V»  [(B®  V)b]  . 

The-  last  term  can  be  wrltteii  in  various  forms: 

r , dB,  6Bx  q^(B,bJ  dB,  dB. 

^ ’ [^5,  ‘ ] = ■sir  ■53T:  = ^ 

^ J ,j.j 

Selecting  the  last  of  these,  sued  substituting  into  (Bl), 


we  find 


it  = hhhj  - hiS?  ^ ‘'m 


^ b2  BB  6B 

„ i _A 

2 ax 


Here  d/dt  = d/dt  + 


T 


j 

V.  Introducing  the  magnetic  stress  tensor 


ij  =n  <h®j  - 5 


(33) 
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f. 


value  of  the  field  at  the  surface  can  be  obtained  by  setting 
X = y = 0 in  (Al)|  the  corresponding  value  (0,0,z)  is  the  field 
at  that  point  on  the  surface  of  the  mantle  which  is  directly  above 
the  center  of  the  upwelling  in  the  core.  'fhis  value  can  be 
expressed  in  the  form 

p3^(0,0,z)  = ■—  7 /%j(p,(p,0)vp2+z2)“V2pdpd^,  (A3) 

o o , . 

Now  p^(0,0j,z)  is  certainly  finite.  hence  there  is  a distance  c 
such  that  “most”  of  the  field  p^(OpOj,z)  comes  from  sources  lying 
rithin  a circle  of  radius  c about  the  origin  in  the  plane  z = 0, 
This  quantity  c could  be  regarded  as  the  "radius  of  the  upwelling," 
so  far  as  its  magnetic  effects  are  concerned.  If  we  now  make  the 
additional  assumption  that  z is  larger  than  c,  we  can  (in  a very 

O 

crude  approximation;  neglect  the  term  p in  the  denominator  of 
(A3).  We  then  obtain 

. c , 

p^(0,05z)~'  / / p.  (pjm.O)pdpd^'^  — ^ I.,  (A4) 

2Hz^  o o 2z‘' 
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Strpplement  * Diffusion-Induction  Equation  for 
the  Energy  Density, 

In  scalar  diffusion  theory,  there  is  often  available  a 
conservation  principle  which  may  be  of  some  use  in  predicting  or 
visualizing  the  course  of  the  diffusion  process.  In  magnetic 
diffusion  theory  there  is  no  field  quantity  which  is  actually 
conserved,  but  one  can  write  the  differential  equation  for  the 
energy  dens  ity  ^^/2  yU  in  the  form  of  a diffusion  equation  with  a 
loss  term. 
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Vi’c  cai:  reduce  this  expression  to  the  slightly  more  symmetrical  form 


2-  4:  = T,  - V . , + 

dt  \ ^ y i.1  ij  ii  ) 


Hi  V 2a 


T as.  3B.  1 

1 i 1 I /X, . \ 

I- 


The  first  two  terms  represent  energy  transfer  from  the  fluid 
motion  to  the  field.  I’he  quantity  V..  is  simply  V • v,  so  bhat 
if  the  fluid  is  incompressible  the  entire  contribution  of  motional 
induction  is  contained  in  the  scalar  product  of  the  magnetic 

stress  tensor  (B3)  and  the  strain  tensor  (B2)  of  the  fluid  motion. 

If  the  fliild  is  at  rest,  the  first  two  terms  on  the  right 
in  (34)  vanich , the  symbol  d/dt  on  the  left  reverts  to  d/dt,  and 

O 

the  equation  reduces  to  a diffusion  equation  for^  /2^  with  a 
loss  term  -(y^/«) (dB^/ox^) (oB^/Bx j ), 

Under  some  circumstances,  the  flux  of through  a surface 
may  be  conserved  during  diffusion  of  the  field.  If  one  can  somehow 
establish,  without  first  solving  the  entire  diffusion  problem,  that 
for  some  specified  curve  G the  right  side  of  the  equation 

It  /^  • ="V'm^ 

5 0 

will  remain  zero  during  a given  time  interval,  then  one  can 
conclude  that  the  flux  through  the  corresponding  surface  S will 
remain  constant  during  that  interval. 
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